Factors influencing mammalian community occupancy in Dhorpatan Hunting Reserve, Nepal

Abstract The decline in mammalian species diversity is increasing worldwide, with areas characterized by high human activities experiencing more prominent effects. Knowledge of spatial distributions of species and factors acting on them is necessary for effective management. We evaluated community‐level occupancy of mammal species in Dhorpatan Hunting Reserve (DHR), Nepal using remote cameras during 15 March–15 June 2022. We used mammal species detections from remote cameras and multispecies hierarchical occupancy modeling to assess the effects of environmental and anthropogenic variables on community‐level occupancy of detected mammal species. We identified a highly heterogeneous mammal species community at DHR with greatest detection probability (0.21) for red fox (Vulpes vulpes) and lowest (0.08) for blue sheep (Pseudois nayaur). Naïve occupancy ranged from 0.31 for giant‐flying squirrel (Petaurista magnificus) to 0.84 for red fox. Mammal community occupancy increased with increasing canopy cover and number of livestock detections, but overall occupancy declined close to human settlements. The findings of this study can be used for developing policy at DHR for the management of mammal species through reducing the potential increase of human settlements or livestock grazing.


| INTRODUC TI ON
Knowledge of species distributions and factors that influence them is a crucial part of ecology and conservation (Gotelli & Colwell, 2001). Recent anthropogenic effects on species' habitats have increased and consequently reduced global biodiversity (Drouilly et al., 2018). These anthropogenic factors impact species occurrences and distributions, which can vary spatially and among species (Karanth et al., 2009;Leweri et al., 2022). Thus, there is a need to monitor and manage species' habitats in response to these threats (Drouilly et al., 2018;Kalle et al., 2013;Leweri et al., 2022).
Community-level information can be gained using multispecies approaches including occupancy modeling, which uses species presence data to gain inferences about their distributions (Balmford et al., 2005;Dorazio & Royle, 2005;Yoccoz et al., 2001). Species-and community-level occupancy is determined by some set of variables (Shmida & Wilson, 1985) which can include anthropogenic factors like settlements, roads, and livestock. For example, mammal community occupancy in China decreased with increasing cattle detections and vehicle traffic but increased with distance to settlement, human detections, and plant productivity (Feng et al., 2021). Similarly in the rangelands of Karoo, South Africa, mammal and ground bird species occupancy increased with increasing livestock detections but decreased with increasing human presence (Drouilly et al., 2018). Environmental factors like canopy cover can also influence occupancy probability of species (Laurance et al., 2008;Whitworth et al., 2019). In Manu Biosphere Reserve, Peru, rainforest mammal species occupancy increased with greater canopy cover and decreased with increasing forest disturbance (Whitworth et al., 2019). Similar effects of canopy cover on species occupancy occurred with nocturnal mammals in African rainforests (Laurance et al., 2008) and mammalian carnivores in Brazil (Regolin et al., 2017).
Under high anthropogenic activities, the probability of encounter between humans and wildlife becomes higher, which leads to conflicts between them (Treves & Karanth, 2003). In these areas, coexistence between wildlife and human should be maintained for their long-term conservation. Assessing animal abundance is prerequisite for their effective management and wildlife conservation (Blanc et al., 2007). For effective assessment, observations of animals should be combined with inferential methods including the potential influence of considered variables (Dorazio & Royle, 2005;MacKenzie & Kendall, 2002;Yoccoz et al., 2001), and these inferences work efficiently only if the system dynamics are generated by manipulative experiments (Fisher et al., 1943;Hurlbert, 1984;Pianka, 1966). Recent advancements in camera-based studies on multiple species have emphasized particular species guilds (Schuette et al., 2013;Stoner et al., 2007), with few efforts to investigate entire communities (Rich et al., 2016), highlighting need to study community-level occupancy. In Nepal, studies of species occupancy have mostly focused on single species (Barber-Meyer et al., 2013;Lamichhane et al., 2021;Sharma et al., 2020;Thapa et al., 2020;Thapa & Kelly, 2017).
Dhorpatan Hunting Reserve (DHR) in Nepal is a highland protected area with human access for natural resource use, which can increase adverse effects between humans and wildlife. However, there have been no investigations of mammalian community response to human activities in DHR. We quantified multispecies occupancy of mammal species and the impacts of habitat variables on their occupancy in DHR. Due to the potential impacts of anthropogenic activities in DHR, we predicted that occupancy of wild mammal species in DHR would be adversely affected by anthropogenic factors such as human settlements and livestock.

| Data collection
For identifying species presence, we deployed remote cameras in Fagune and Barse blocks of DHR during 15 March-15 June 2022. We first created a grid of 290 plots 1 km 2 cells for camera deployment.
We deployed a single camera (Stealth Cam STC-G45NG) in each of 34 randomly selected cells (116 km 2 overall area), with a spacing of about 1 km between cameras in adjacent cells. We placed cameras 50-100 cm above ground and typically near trails, water bodies, or other areas where species were likely to be detected. We set cameras to take three images for each detection with a 30-s delay. This study took each 14 days survey period as a single sampling occasion across the study.
We measured canopy cover, distance to nearest human settlement, and number of times livestock were detected at each camera site. We recorded forest canopy cover using Gap Light Analysis Mobile Application (GLAMA; Tichý, 2016) in each camera trap location. We established a 10 × 10 m plot at the center of each camera trap location and measured the canopy cover from four corners and the center of each plot. We measured distance to nearest human settlement using a measuring tape if less than 200 m and Google Earth when the distance exceeded 200 m. We determined the number of times livestock were detected using images from each camera. Because large herds generally occur in DHR, we counted sequential detections (at 30-s intervals) as a single detection when it was apparent a single herd was passing the camera. We classi-

| Data analysis
We filtered, arranged, and summarized the data using the package tidyverse (Wickham & Wickham, 2017) in R Program (R Core Team, 2022), then created a detection matrix with the number of detections for each species across occasions. Before analysis, we performed the correlation and collinearity analysis to identify whether the independent variables were highly related or not. The analyses revealed none of the selected variables were correlated (Pearson's correlation |r| < .7) (Dormann et al., 2013)

| Hierarchical multispecies occupancy
We conducted the hierarchical analysis following Royle and Dorazio (2008) in which the Beta was replaced prior for specieslevel ψ k : and the community-level hyperpriors were added as We subjected the analysis to Markov Chain Monte Carlo (MCMC) simulations to obtain the posterior distributions. We ran 50,000 iterations for each of three chains then visually inspected whether the chains mixed well. We confirmed mixing using trace plot diagnostics and Rhat values. If the optimal value of Rhat was <1.1, we did a further step in analysis by adjusting the number of iterations until the Rhat became <1.05. After the successful convergence of markov chains for all parameters, the value of detection probability (p) and naïve occupancy (ψ) including mean and standard deviation (SD) of p and ψ were obtained.

| Effect of covariates
As occupancy probability varies across sites, we calculated occupancy for each species k at site i using the following formula: where the coefficients (β) differ among species. We used species level priors as: where μ x and σ x are the coefficients of random variable drawn from a normal distribution with mean and SD to be estimated. Unlike previous models, we incorporated correlation between ψ and ρ, where the intercept β 0,k , is the probability of occupancy of species k at a site with a given combinations of variables (Devarajan et al., 2020).
For community-level hyperpriors, we used a uniform Beta (1,1) prior to the probability and converted it to the logit scale where: ∼ Beta(1, 1) and 0 = logit 0 . We used a uniform prior for the SD as: Σ 0 ~ Uniform (0,5). The coefficient value at logit scale will be around ±5, therefore, μ x = Uniform (−5, 5) and the SD is σ x = Uniform (0,5). For hierarchical analysis, we used an adaptive MCMC with 50,000 iterations, three chains, 1000 adaptations, and a burn in of 1000. We performed occupancy analyses using Just Another Gibbs Sampler (Plummer, 2003) and R Program (R Core Team, 2022) using packages coda (Plummer et al., 2006) and jagsUI (Kellner et al., 2019).

| RE SULTS
We detected 15 mammal species (83 total detections) from eight families and four orders in 1530 camera days. One species (red panda) is categorized as endangered, four species as vulnerable, one near threatened, and nine as least concern in the IUCN Red List

| Effects of covariates
Mammal community occupancy appeared to decrease when closer to human settlements (1.451 ± 1.107) but the effect was not significant.
The occupancy of all species increased with increasing livestock detections except for Himalayan goral (intercept = 1.70 ± 2.95) which remained relatively constant (Figure 4).
The greatest increase in occupancy with number of livestock detections was observed for red panda (intercept = 4.14 ± 2.50).
We observed a positive impact of livestock detection on occu- Yellow-throated marten (intercept = 0.01 ± 1.85) was observed as distance from settlements increased.

| DISCUSS ION
The study confirmed the heterogenous mammal species composition in the study area and anthropogenic effects such as livestock presence on their habitat is rampant. Basically, we found overall avoidance of humans (i.e., settlements), with mammals selecting for greater forest cover. It means the anthropogenic activities are altering (largely negatively) mammal communities. The high value of SD for the covariates means that the study plots also were highly heterogeneous in the variable composition. We also identified a lesser to no effect of the covariates on the species like red fox, barking deer, and Himalayan goral. These species are widespread species in Nepal and are more generalist in their habitat selection (Jnawali et al., 2011); therefore, no major impact of variables was observed. Number of livestock detected and canopy cover were observed to exert a positive impact on mammalian community occupancy in our study human settlement exerted a negative impact.
Number of livestock detection was observed to influence the most on species occupancy positively with a sharp increase in response.
It might be due to common sharing of grazing grounds between livestock and wild herbivores or presence of both prey and carnivore species in the same habitat (Kalle et al., 2013). The DHR is one of the pastureland; therefore, local people or people near to the reserve leave their livestock for grazing which might support the coexistence of wild mammal species and livestock. The highest positive impact of livestock detection was observed on red panda which might be due to sharing of bamboo as food by red panda and livestock (Sharma et al., 2014). The positive impact of herbivore on red panda presence was also reported in protected area system of Nepal (Acharya et al., 2018). The increased occupancy of wild boar with near to F I G U R E 2 (a) Detection probability (p) and (b) Naïve occupancy (ψ) for each of the 15 mammal species across five replicate surveys in Dhorpotan Hunting Reserve, Nepal, 2022.

TA B L E 2
Community-level summaries of the hyperparameters for occupancy hypothesized to influence the occupancy probabilities of mammalian community in Dhorpotan Hunting Reserve, Nepal, 2022.  human settlement in this study indicates the potentiality of this species for human-wildlife conflict in DHR as found in other protected areas (Pandey et al., 2016). However, in case of Himalayan goral, the relatively less impact of livestock presence might be it being a widespread species in the study area. The relatively low detection rate of blue sheep within the study area could potentially be attributed to a decline in population size in recent years, as well as their known preference for terrain slopes as a habitat (DHR, 2019).
The increase in forest canopy cover supports to increase in the occupancy of all species. It might be due to their preferences toward dense habitat potentially to avoid human disturbances and hiding place with predators (Laurance et al., 2008;Regolin et al., 2017;Whitworth et al., 2019). The dense vegetation provides a better hideout for prey species as well as ambushing spots for predators (Monroy-Vilchis et al., 2009). Himalayan serow was the species most positively impacted by the percentage of canopy cover. The preference of dense vegetation covers by Himalayan serow was also observed in Kanchendzonga Biosphere Reserve, India (Bhattacharya et al., 2012), which might be due to it being prey species and preferring dense vegetation to avoid predators.
Occupancy of all species increased with the increase in distance from human settlement except for Yellow-throated marten. The increase in occupancy of all species with the increase in distance from human settlement generally relates to the increase of anthropogenic disturbance as human settlements offer different threats to wildlife like poaching as well as guard dogs (Cavada et al., 2019;Salvatori et al., 2022;Schuette et al., 2013). The decrease in occupancy of species with proximity to the settlement due to the hunting and poaching pressure is rather common observation in case of ungulates (Soh et al., 2014). Carnivores also prefer to roam away from settlements and are observed to avoid possible human encounters to avoid poachers as well as retaliation from local farmers (Drouilly et al., 2018;Kalle et al., 2013;Pia et al., 2013). The decrease in occupancy of Yellow throated marten with increase in distance from human settlement correlates to its prey preference as it prefers small livestock prey species like avian livestock (Baral et al., 2021) and human associated small mammal prey species like rats (Rattus rattus; Parr & Duckworth, 2007;Zhou et al., 2008).
The occupancy of barking deer, red fox, and Himalayan goral remained relatively stable across all variables, which might be due to higher species abundance as well as the ecological aspects of the species in the study area. These species are observed to have widespread distribution range across Nepal (Jnawali et al., 2011). Higher site occupancy for Barking deer was observed in the area where low impacts of variables are noticed (Letro et al., 2022). The occupancy of red fox was also high across different habitats in Tieqiaoshan Nature Reserve, China, with less extinction probabilities in sites (Vitekere et al., 2020) as well as relatively stable impacts of associated variables (MacDougall & Sander, 2022). Himalayan goral is also one of the species with higher occupancy across the study area and was observed to be less impacted by the variables involved. This might be due to the lesser preference of elevation by the species as well as lesser selectivity of habitat by the species (Bhattacharya et al., 2012).

ACK N OWLED G M ENTS
We thank University Grants Commission for financial support. We thank Department of National Parks and Wildlife Conservation and DHR for providing research permission.

F I G U R E 4
Effects of covariates number of livestock detected, canopy cover, and distance from human settlement on the occupancy of 15 mammal species, Dhorpatan Hunting Reserve, Nepal, 2022.

CO N FLI C T O F I NTE R E S T S TATE M E NT
Authors declare no conflict of interest.

FU N D I N G I N FO R M ATI O N
University Grants Commission, Nepal (77/78-S&T-38).

DATA AVA I L A B I L I T Y S TAT E M E N T
All the relevant data used in this study will be archived in Dryad https://doi.org/10.5061/dryad.6m905 qg3s.